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2 ! Abstract 

k>Y We systematically study K + observables in nucleus-nucleus collisions at 1-2 A GeV within the 

Boltzmann-Uehling-Uhlenbeck (BUU) transport model. We compare our calculations with the 
KaoS data on the kaon multiplicities and spectra. In addition, the kaon collective flow is computed 
and compared with the FOPI and KaoS data. We show, that the elliptic kaon flow measured 
■ recently by the KaoS Collaboration is best described by using the Brown-Rho parametrization of 
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I. INTRODUCTION 



Since about 20 years strangeness production in heavy-ion collisions is a hot topic of 
theoretical and experimental studies. Due to rather high energy thresholds in NN collisions 
(E beam = 1.58 GeV for NN -> KAN and E beam = 2.5 GeV for NN -> NNKK) the 
secondary processes AN — > KYN, tcN — > ITK and 7rV — > A'iV - which require high baryon 
density - are important in the case of nucleus-nucleus collisions at 1-2 A GeV. Moreover, 
due to the relatively low (~ 10 mb) KN scattering cross section and the absence of the 
absorption channel of a kaon on a nucleon in strong interactions, nuclear matter is practically 
transparent for kaons 1]. Thus, the kaon yield is a good probe for the nuclear equation- 
of-state (EOS). This idea has originally been proposed and tested in BUU calculations by 



Aichelin and Ko in Ref. [2j. Recently it has been shown by Fuchs et al. |3] within the 
Tubingen QMD model, that the ratio of the kaon yields in Au+Au and C+C collisions 
plotted vs the beam energy favours a soft EOS (K=200 MeV) when comparing with the 
KaoS data J|. This ratio seems to be weakly sensitive to the experimentally not measurable 
AN — > KYN cross section and to the choice of the kaon potential and serves as a good 
probe for the nuclear EOS. 

Another direction of studies is related to the kaon and antikaon propagation in the nuclear 
.edien, EflflflflHHQ.Fo, thrae studies the K and K m ean fl e!d potent p,a y a 
crucial role. In the lowest order approximation to the chiral Lagrangian the kaon (antikaon') 
potential has an attractive scalar and a repulsive (attractive) vector part jfil IrIIi^I]. This 
leads to a weakly repulsive potential {Uk{po) — 7 MeV, where p — 0-17 fm -3 ) for kaons 
and a strongly attractive potential (^(po) — —100 MeV) for antikaons. The kaon P, |(| 
and the antikaon in-plane flows and the kaon azimuthal distributions [?| are strongly 
influenced by the K and K potentials. The authors of Ref. I| have shown within the HSD 
model, that a repulsive potential (Uk(po) — 30 MeV) seems to be needed for description of 
the first FOPI data UM on the K+ in-plane flow. In the most recent HSD analysis of Ref. 

n 

|12| . where different kaon potentials were tested, the one given by the chiral perturbation 
theory in the relativistic Hartree approximation (Uk{po) — 20 MeV) was found to give the 
best agreement with the kaon flow data. The same value of the kaon potential at normal 
nuclear matter density (Uk(po) = 20 ± 5 MeV) was reported in a recent CBUU analysis of 
the K + production in proton- nucleus reactions by Rudy et al. 
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In the approaches of Refs. 



BBBBQQ 



the spatial components of the vector field in 



the kaon potential were neglected which should be approximately valid in a central zone of 
a symmetric colliding system where the barvon current disappears. The spatial components 
were taken into account by Fuchs et al. |l7( within the Tubingen QMD model and led to a 
much smaller negative flow of K +, s. This is a consequence of the cancellation effect of the 
repulsive time component of the kaon vector field by the Lorentz force generated by the space 
components of the kaon vector field However, the calculations of Ref. [loj . where the 
spatial components have also been taken into account, have produced a larger negative kaon 
low, which disagrees with Ref. The latest calculations by the Tubingen QMD group 

11 1 have corroborated their earlier analysis and demonstrated that the new FOPI data 
on the kaon in-plane flow are best described by using the kaon potential given by the 
Brown-Rho (BR) parametrization (Uk(po) — 30 MeV). At p < p the kaon potential in 



the BR parametrization is close to the one in the impulse approximation(c.f. Refs. |6l l20|). 
Recent self-consistent calculations of Refs. J21I Q show even stronger repulsion for kaons 
(U K {po) = 36 MeV in Ref. Q and U K (po) = 39 MeV in Ref. Q). 

The K ± azimuthal distributions at midrapidity (squeeze-out) have also been studied 
within an earlier version of the Tubingen QMD model in Ref. . It has been concluded in 
23I ] , that the K + squeeze-out is caused mainly by the repulsive K + potential. However, the 
quantitative agreement with experimental data for Au+Au at 1 A GeV 24\ was achieved 
only with the static kaon potential in [23]. The Lorentz force, as has been concluded in j^j], 
destroys the agreement with the data by reducing the K + squeeze-out signal strongly (we 
will come back to this point later on). 

The present work is an attempt to describe the K + data on mult i plic ities and phase 
space distributions at the beam energies of 1-2 A GeV J L 



the basis of a BUU model 



30 
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3J. The BUU model 



IJ18J24J25I, 

hd, ki M i 



26 



27 



28|, |29J on 



includes a large set 



of the baryonic resonances (see the next Section). Due to the channels NR — > KYN and 
NR — > KKNN, where R stands for a nonstrange baryon resonance, the resonances enhance 
strangeness production. Thus, despite of many previous transport theoretical studies on the 

BBBBBBBQHI 



12, 



33|,|3J) it 



strangeness production at SIS energies (c.f. Refs. 

would be interest ing to confront also our calculations with experimental data, since, e.g., the 



QMD models 



esting 

m 



JU1IJ33L34, 



and the HSD model 



351I36J propagate only A(1232) and A"* (1440) resonances 



propagates only A(1232), A*(1440) and A*(1535) resonances 
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(see also Ref. |38[ for the comparison of the different transport calculations). 

The main purpose of our study is to clarify whether a kaon potential is actually needed 
to describe the data and, if so, how strong it must be. In particular, we will analyse the 
recent KaoS data on the kaon azimuthal distributions. 

In Sect. II a brief description of the BUU model is given. Sect. Ill contains the results of 
numerical calculations. In Sect. IV we summarize our results and draw some conclusions. 



II. THE BUU MODEL 

Our calculations are based on the BUU model in the version described in Refs. 



ail 0, Q- 

The model explicitly propagates all N* and A resonances that are rated 
with at least 2 stars in the analysis of Ref. j^, which includes the N* states Pn(1440), 
D 13 (1520), Sn(1535), 5 n (1650), D 15 (1675), Pi 5 (1680), P 13 (1879), F 17 (1990), Gi 7 (2190) and 
the A states P 33 (1232), P 33 (1600), S 31 (1620), 5 3 i(1900), P 35 (1905), P 3 i(1910), P 35 (1930), 
P 37 (1950), P 35 (23 50). Also the S = -1 baryons Y = A(1116), S(1189) and Y* = 5 i(1405), 
P 03 (1520), P i(1600), 5 01 (1670), A> 3 (1690), S i(1800), P i(1810), P 05 (1820), P 05 (1830), 
P 03 (1890), G 07 (2 1 00), P 05 (2110), P 13 (1385), P u (1660), D 13 (1670), S n (1750), #15(1775), 
Pi5(1915), Pi 7 (2030) are propagated explicitly in the model. The y*-resonances are the 
intermediate states of the strangeness exchange reactions tcY <-> Y* <-> KN . 

In the meson sector, the following particles are propagated: ir, r], p, a, lu, rj', 0, rj c , J /if), 
K, K, K* , K* . Also the cascades, charmed baryons and mesons are propagated, which are, 
however, irrelevant degrees of freedom at SIS energies. 

The m ode. has bee, s— y applied to 7 and e induced .eaet.ons on nnclei HQ 
and to the description of the collective nucleon flows in heavy-ion collisions at 0.1-2 A 
GeV j41(. The reproduction of the pion abundancies in heavy-ion collisions at SIS energies 
requires, however, to apply the in-medium reduced resonance production/absorption cross 



sections NN <-> NR 
be found in Refs 



32 



42 



_ ._ . We will drop the detailed description of the model which can 

rrnm 

• 30J, |3JJ, |32J concentrating here only on some novel features related to 



strangeness production. 

For kaon production in baryon-baryon collisions we include the following channels: NN — > 
NYK, NN -> AYK, AN -> NYK, AN -> AYK and AA -> AYK, where A = P 33 (1232) 
and Y = A(1115) or £(1189), with all possible isospin combinations of the incoming and 
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the outgoing particles. The channel AA — > NYK is not included, since it is possible only 
through an exchange by an on-shell pion j^, which is already included in BUU via the 
consecutive A <-> Ntt and tiB — > KY processes, where B = N or R. The isospin-dependent 
parametrizations for the cross sections of kaon production in the NN, NA and AA collisions 
are taken from Ref. |43j]. In addition, we allow to produce a kaon-hyperon pair in a collision 
between any two nonstrange baryons: cross sections of these processes are obtained from 



the cross sections of Ref. 



43j by replacing an incoming N* resonance by the nucleon and an 



incoming higher A resonance by the A(1232) for the same ^/s. The KK pair production in 
a baryon-baryon collision BB — > NNKK is also included via cross sections parametrized 
in Ref. Q 

The pion-baryon collisions ttB — > YK also contribute quite significantly to the 
strangeness production in heavy-ion collisions. If the incoming baryon B is a nucleon or 
a A(1232)-resonance, the cross sections from Ref. 45| are applied. If B = N* or higher A, 
it is substituted by the nucleon or A ( 1 2 32 1 1 -re sonance , depending on the isospin, and then 
the corresponding cross sections from Ref. [45] are used. The ttB NKK process is taken 
into account with a cross section parametrized in Ref. j^fj. The kaon elastic, including 
charge exchange, and inelastic scattering processes KN — > KN and KN —>■ KNir are also 
taken into account: their cross sections are fitted to the data (see Ref. j^). 

Since in the present work we concentrate on the kaon production at energies below 2 A 
GeV, the kaon production channels BB — > NYK and ttB — > YK are those relevant for our 
study. The reactions including antikaons play, practically, no role for kaon production (see 
Tabledbelow). The kaon production is treated perturbatively and the FRITIOF mechanism 
is switched off in this work, since the kaon multiplicity per nucleus-nucleus collision is still 
quite small at the considered beam energies (see Fig. El below) . 

In our calculations the nucleons are prop agated in a Skyrme-like mean field including the 
momentum-dependent part (c.f. Ref. 4l|). Most of the calculations are done with a soft 
momentum-dependent mean field (SM, K = 220 MeV). When, for comparison, also a hard 
momentum-dependent mean field (HM, K = 380 MeV) is used, this will be mentioned in 
the text below. It is assumed that all the nonstrange baryonic resonances experience the 
same mean field as nucleons. For the hyperons, according to the fraction of the nonstrange 
quarks, we apply the nucleon mean field multiplied by 2/3 P. I37I l^. 
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The K single-particle energies are expressed as 



c4(k) = ±V" + ^k* 2 + m*£ , (1) 

k* = k =F V is the kaon kinetic momentum, V 1 = (V°, V) is the kaon vector field and m* K 
is the kaon effective (Dirac) mass. According to Refs. |llL Il9j. the kaon effective mass and 
the kaon vector field are given by the following expressions : 



m' K = Jm' x - ^fp, + V- , (2) 

V = • (3) 

where mx = 0.496 GeV is the bare kaon mass, p s and j M are the baryon scalar density and 
the four-current, respectively. The parameters which appear in Eqs. are ^kn = 0.450 

GeV — the kaon-nucleon sigma term, f n = 0.093 GeV — the vacuum pion decay constant 
and ft = y/(L6f n — the in-medium pion decay constant at normal nuclear matter density 

Using the in-medium pion decay constant rather than the vacuum one in Eq. leads 
to the desirable consequence fl3 | . that the resulting kaon vector potential Q is just 1/3 of 
the nucleon vector potential given by the relativistic mean field model. In the scalar term 
oc p s (Eq. 0) the vacuum pion decay constant is used, since for kaons the higher order 
(range) term in the chiral expansion cancels an effect of the in-medium pion dec ay c onstant 

191. Thus, 



in this case . The cancellation, however, does not take place for the antikaons 
Eq.(j2J) can be considered for the antikaons as a phenomenological parametrization, which, 
nevertheless, leads to reasonable values of the K potential (see below). 
The kaon (antikaon) potential U K r^\ is defined as 



U K{R) ^)=^)-^ + ml. (4) 

Following Ref . ^j] , we will denote the kaon (antikaon) potential with parameters determined 
above as the BR potential which will be used as the default potential in our calculations. 
However, as in 11], we will perform for a comparison also some calculations applying the 
kaon (antikaon) potential with T^kn = 0.350 GeV and with the free pion decay constant f n 
instead of /* in the vector potential © which we call a Ko-Li (KL) potential 0. 

Both choices of the kaon potential, BR and KL, are shown for the static case (k = 0) in 
the upper panel of Fig. ^ as functions of the baryon density. The BR potential is strongly 



repulsive (Uk(po) = 32 MeV), while the KL potential is much weaker (Uk(po) — 6 MeV). 



We have to point out here, that our potentia 
than the corresponding potentials from Ref. 



s are slightly more repulsive at higher density 



111 ] due to different versions of the relativistic 



mean field model used to evaluate the scalar density. For an orientation, we show in the 
lower panel of Fig. [T] the scalar density as a function of the baryon density which is given 
by the NL2 model (49| applied in our calculations. In the case of antikaons, the BR (KL) 
parametrization produces £/g-(po) = —144 MeV (-100 MeV). 

In the BUU implementation we solve the Hamiltonian equations of motion for the kaon 
(antikaon) test particles, where the Hamilton function is the single-particle energy (0) which 
now implicitly depends also on space and time via the vector field V^ir, t) and the effective 
mass m* K (r, t) : 

&4(k,r,t) k* 



<9k E* ' 

du%{k,v,t) m* K dm* K dV k* a dV c 



(5) 



dr E* dr dr E* dr 



where E* = wk* 2 + m*^. These equations of motions are completely equivalent to those 



in the covariant form derived in Refs. 
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\v\ . The last term in the r.h.s. of Eq.(JHJ) is 

a velocity-dependent (Lorentz-like) force caused by the spatial components V of the kaon 
vector field. 

Potentials shift the particleproduction thresholds in nuclear medium. We take this effect 
into account, following Ref. |3j, by replacing the bare center-of-mass (cm.) energy y/Sf ree 
in the argument of the cross sections crBB^ s trangeness(v^7^) and ^^strangeness (v/s/^) by 
a corrected quantity. E.g., for the process iVW — > NYK we replace ^/Sf ree by \fs — fh^ — 
rhy — fhic + rriN + my + tuk, where is the total in-medium cm. energy of colliding 
particles including their mean field potentials, rhx = £x(Px m ' = 0) (X = N, Y, K) are 
the in-medium masses of the particles defined as their energies at rest in the cm. frame of 
colliding particles jscj. 

In the simulation of a three-body phase space for the outgoing particles BYK (B = N 
or A (1232V here) we use the kaon momentum distribution in the BYK cm. frame proposed 
in Refs. [jj, kl| : 

dW K ~ (— Vf 1 -— V . P) 
where p max = [{s — (tub + fhy) 2 + fh 2 K ) 2 /As — fn 2 K ] l l 2 is the maximal kaon momentum in the 



BYK cm. 
form 
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rame. The angular distribution of the produced kaon is chosen in the empirical 
52( : da/dcos Q c .m. oc (1 + a cos 2 C . TO .) with a = 1.2 5^|. These two modifications 
lead to a softer kaon pi ab spectrum with respect to the simulation using an ideal Dalitz 
3-body decay 

As it was demonstrated in Ref. |32j, taking into account the Dirac masses of the baryons 
reduces the cross sections NN «-> NR in nuclear matter strongly, which brings our BUU 
calculations in a better agreement with the experimental data on pion multiplicities. In the 
calculations of the present work, if opposite is not stated explicitly, we use the in-medium 
cross sections NN <-> NR and NN «-> NNn 0] computed with the Dirac masses from 



the NL2 model |49|. The exchange pion collectivity effect and vertex corrections in the 
NN <-» AA(1232) cross sections included in the calculations of Ref. 3^ are neglected in 
the present work for simplicity. For brevity, we will call below a calculation employing the 
SM nucleon mean field, the BR kaon potential and the in-medium iViV «-> NR cross sections 
the standard one. 



III. NUMERICAL RESULTS 
A. Kaon production channels 

Before starting the comparison with experimental data, we have looked at the time evolu- 
tion of the kaon number produced by the different channels, which is displayed in Fig.|2k,b,c 
for the central collisions Au+Au at 0.96 A GeV, Au+Au at 1.48 A GeV and Ni+Ni at 1.93 A 
GeV, respectively. Contributions of the channels to the total kaon number are summarized 
in Table [fl In Fig. |21 and in Table E U B" in the initial state denotes N or R, while in 
the final state denotes N or A(1232). The main kaon production channels are RN — > BYK, 
NN — > BYK and irB — > KY. Other channels contribute all together on the level of ~ 10%. 

At the lowest energy of 0.96 A GeV the direct channel NN — > NYK is deeply subthresh- 
old and, thus, its contribution to the total kaon number is small. With increasing beam 
energy the contribution of the NN — > NYK channel grows quickly, so that at 1.93 A GeV 
this channel becomes already the main one. The RN — > BYK channel dominates at 0.96 
and 1.48 A GeV. It is also the second important at 1.93 A GeV. 

The kB — > KY channel is the second important at 0.96 A GeV, but its contribution 
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decreases with increasing energy. This is quite natural, since a channel caused by secondary 
particles is energetically favourable at subthreshold energies. Another reason is that the 
pions are produced mainly in resonance decays R — > Nn. The time scale of a nucleus-nucleus 
collision gets shorter with increasing beam energy. Therefore, more and more resonances 
will decay at the final low baryon density stage, when the probability of the pion-baryon 
collisions is small. 

By comparing the time evolution of the kaon production (Fig. |2k,b,c) and of the central 
baryon density (Fig. EH) we see that the largest rate of the kaon production takes place at 
about the time of the maximum compression. Thus, the kaon production delivers a signal 
from the highest compression stage of a heavy-ion collision without much distortion from 
the later stage, since kaons are not absorbed in the nuclear medium. This important kaon 
property has given, in particular, an opportunity to determine the nuclear EOS from the 
kaon yields |2j, |3j . 



B. Kaon multiplicities and spectra 

First, we check whether our model is able to reproduce the pion and kaon total mul- 
tiplicities. Fig. El shows the inclusive multiplicities of 7r's (n = tt~ + 7r° + tt + ) and K + 's 
per projectile nucleon for the systems Au+Au and C+C as functions of the beam energy 
in comparison to the data from Ref. j^J. In the calculations the particle multiplicities were 
impact parameter weighted in the region b < 14 fm for the Au+Au collisions and in the 
region b < 5 fm for the C+C collisions which corresponds to the geometrical cross sections. 
The pion multiplicity for the C+C system is well reproduced, except for the points at 1.8 
and 2.0 A GeV, where we underpredict the data due to neglecting the string (FRITIOF) 
mechanism of particle production. In the case of Au+Au collisions our calculation produces 
too many pions at lower beam energies in spite of using the in-medium reduced NN <-» NR 
cross sections. Applying faster dropping Dirac masses of the baryons with nuclear density, 
e.g. given by the NL1 model as well as taking into account the pion collectivity and 
vertex corrections in the iViV <-> iVA(1232) cross sections [3^ would lead to a better de- 
scription of the pion multiplicity in the Au+Au case, however, at the cost of too low pion 
multiplicity for C+C collisions. The K + multiplicities are rather well reproduced for both 
systems Au+Au and C+C in the full SIS energy range. 



9 



Fig. |3] shows the laboratory momentum kaon spectra for the Au+Au collisions at 0.96 A 
GeV in comparison to the data [4]. The calculation without kaon potential (Fig. E^, dashed 
lines with open up triangles) overestimates the kaon production strongly. The calculations, 
which include the kaon potential (Fig. E^, BR — solid lines with open circles, KL — dotted 
lines with open squares) lead to a reduced lower momentum part of the spectra in better 
agreement with the data. The high momentum tail of the K + spectrum still remains too 
high in the calculations with the kaon potential. 

In Fig. Et> we compare our standard calculation with the calculation employing the vac- 
uum NN <-> NR cross sections. We see that dropping the in-medium corrections to the 
iViV <-> NR cross sections results in about 50% larger kaon production cross section. 

In order to see the origin of the produced kaons, we performed in Fig.Efc a channel decom- 
position of the laboratory momentum spectrum 5^ . The largest contribution is provided 
by the pion-baryon collisions {ttB). The nucleon-A(1232) and nucleon-higher resonance 
channels have also big contributions, comparable to the channel. This explains also a 
sensitivity of the kaon production cross section to the choice of the iViV <-> NR cross sec- 
tions (see Fig. Eb). The nucleon-nucleon collisions contribute only about 10% to the total 
kaon production cross section, since, at = 0.96 A GeV, they are mostly below the kaon 
production threshold. 

In Fig. EJi we present a calculation with the HM nucleon mean field, which produces less 
K + 's (dashed line with open diamonds) than the standard calculation employing the SM 
mean field (solid line with open circles). The pion off-shellness effect j3| (dotted line with 
open squares) does not increase the K + multiplicity strong enough to bring the calculation 
with the HM mean field to agreement with the data. 

Fig. El shows the laboratory momentum kaon spectra for the system C+C at 2 A GeV 
in comparison with the data from Ref. 5^]. It is interesting that for this light system we 
reproduce rather well the slopes of the experimental spectra. A similar result was obtained 
earlier in Ref. ^j] in a calculation with an anisotropic angular distribution of the produced 
kaon in a NN collision. 

Fig.lHlshows the cm. kinetic energy spectra of K +, s from Au+Au collisions at 1.5 A GeV 
for various event centrality classes. For the most central collisions, the calculation without 
potential (upper dashed line with open triangles) clearly overestimates kaon yield and also 
has a too steep slope. Including kaon potentials (BR — solid lines with open circles, KL - 
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dotted lines with open squares) reduces the kaon yield and reduces the steepness of the slope 
inproving, thus, an agreement with the data for central events. With decreasing centrality 
the slopes of the spectra calculated with and without potentials get similar. This is expected 
since in peripheral events the compression is less and, therefore, any influence of the kaon 
potential is reduced. For the most central events we observe that the calculation employing 
the BR potential describes the data best, while for the peripheral events the calculation 
without potential is in best agreement with the data. 

The qualitative trends presented in Fig. can be better visible if one fits the spectra as 
follows [28]: Ed 3 a /dp 3 = C ■ E ■ exp(—E/T), where E = Ef.™ + m K is the total energy 
of a kaon in the cm. system. In the upper panel of Fig. [7| we show the inverse slope 
parameter T of the K + cm. kinetic energy spectra as a function of a participant number 
Apart for the Au+Au system at 1.5 A GeV. The participant number was determined for each 
impact parameter from the geometrical overlap of colliding nuclei assuming sharp nuclear 
surfaces. As we already saw in Fig. |H1 the difference between the inverse slope parameters 
from the calculations with and without kaon potentals is small at the peripheral collisions 
and increases with the collision centrality. The calculation with the KL potential provides 
the best description of the inverse slope parameter for all A part . 

The lower panel of Fig.[7|shows the K + multiplicity per participating nucleon as a function 



of A part . The calculation without potential is above the data by a 
potential is not strong enough to get the correct multiplicities |5J 



actor of two. The KL 
. Only using the BR 



potential reduces the kaon multiplicity to a good agreement with the data. 



Figs. |H1 and El show, respectively, the K + transverse mass (mj = \Jpx + m 2 ) and cm. 
rapidity spectra for Ni+Ni collisions at 1.93 A GeV in comparison to the data Q. As 
in the case of Au+Au collisions at 1.5 A GeV we observe that the calculation without 
kaon potential overpredicts the kaon yields and produces too steep slopes at low m T 's for 
the central collisions (left panels of Figs. EI and The BR potential provides the most 
reasonable description of the data for the central collisions. However, the steepness of the 
slopes is somewhat underestimated by the calculation with the BR potential. The slopes 
are, again, best described by the calculation with the KL potential. In the case of peripheral 
collisions (right panels of Figs. |H1 and OJ) we observe that the calculations without potential 
and with the KL potential produce an equally good description of the data, while the BR 
potential leads to a slight underestimation of the kaon yield. 
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C. Collective flows 



Extraction of collective flows requires a knowledge of a reaction plane, i.e. the plane 
which is parallel to the impact parameter and to the beam momentum j^q . In theoretical 
calculations the reaction plane is obviously known. Experimentally, however, the reaction 
plane is usually taken parallel to the difference of transverse momenta of charged particles 
- which are mostly protons at SIS energies — in the forward and backward hemispheres 
in the cm. frame of colliding nuclei. Thus, a collective flow carries an information about 
the correlations between the particle under study and the remaining protons. This makes 
the collective flow observables extremely useful to constrain the mean field potentials of the 
particles (c.f. Refs. QflflflQ for K and K flows). 

Fig. 1 101 shows a mean transverse momentum projected on the reaction plane as a function 
of a normalized rapidity = (y /y pro j) c . m . for A hyperons (upper panel) and for protons 



(lower pane 



'or the system Ni+Ni at 1.93 A GeV. For comparison, we have selected the 
FOPI data 2Jj|. The proton flow (= d(p x )/dY^ ' at = 0) is quite well described. In the 
calculation of the A flow we have taken into account both A's and S°'s, since in the data 2^ 
both kinds of the A hyperons, those primary and those originating from S° — > A + 7 decays, 
are indistinguishable. For simplicity, we have neglected, however, a difference between the 
momenta of the decaying S° and the outgoing A. A correction for a recoil momentum due to 
the photon emission should slightly reduce the calculated flow. Nevertheless, we observe a 
rather good agreement with the data on the A-flow. This supports our choice of the hyperon 
mean field (see Sect. II). 

It was shown in Ref. Q| that the FOPI data ^| on the K + in-plane flow can be reason- 
ably well described using the BR parametrization of the kaon potential. Before discussing 
an out-of-plane flow, we will also compare our calculations with the FOPI data [tsI. I27I on 
the K + in-plane flow. 

Fig, im shows the K + in-plane flow for Ni+Ni collisions at 1.93 A GeV in comparison to 
the data from Ref. . In full agreement with earlier calculations of Ref. Q| we observe 
that the BR parametrization describes the data best (the solid line with open circles). The 
KL parametrization (dotted line with open squares) does not give an enough repulsion to 
get the negative flow. Neglecting the space component V of the kaon vector field in both 
calculations, i.e. using static potentials, results in too negative flow for both the BR (solid 



12 



line with open pentagons) and the KL (dotted line with open diamonds) parametrizations. 



This was also pointed out earlier in The compensation of the repulsive static potential 
by the V field is so strong, that, e.g. the KL parametrization gives practically the same 
(positive) flow as the calculation without kaon potential (dashed line with open triangles) ! 

An azimuthal distribution with respect to the reaction plane can be represented by a 
Fourier expansion : 

dN 00 

— (0) oc 1 + ]T 2v n cos{n<f)) . (8) 

a< P n=l 

The first two coefficients, V\ and v 2 , in Eq.fjHJ) are called directed and elliptic flow, respec- 
tively, and are expressed as follows : 

V! = (cos0) = (p x /pt) (9) 
v 2 = (co S (2<f>)) = {(p 2 x -p 2 y )/p 2 t ) . (10) 

By neglecting the terms with n > 4 in (jHJ, the elliptic flow can be related to a ratio of 
particle numbers emitted out (0 = ±90°) and in (0 = 0° and 180°) the reaction plane : 

dN/d(j)(90 ) + dN/d(j)(-90 ) _ l-2v 2 

dN/d(f)(0 ) + dN/d(f)(180 ) ~ 1 + 2v 2 ' ^ ' 

Fig. fT21 shows the directed flow V\ as a function of the transverse momentum p t for protons 
around the target rapidity from semicentral Ni+Ni collisions at 1.93 A GeV. We, first, 
considered all the protons including those bound in the target spectator remnant (dashed 
line with open circles). Then we selected only the protons separated from other particles 
by a critical distance d c > 3 fm (solid line with open squares. The proton directed flow is 
negative in agreement with Fig. El We observe a sensitivity of V\ at low p t 's to the selection 
of protons: the calculation including only the separated protons gives a somewhat smaller 
absolute value of V\ than the calculation including all protons. The directed flow at large 
Pt's is independent of the proton selection procedure, and we fail to describe the data here. 
We ascribe this problem to a too hard momentum dependence of the nucleon mean field at 
large momenta 41], which pushes the high-momentum protons too early from the system, 
before the directed flow develops. We do not expect, however, that this drawback influences 
kaon propagation, since the high-p 4 nucleons are not abundant and their contribution to the 
kaon mean field is small. 

Fig. EH displays the directed flow of kaons near target rapidity as a function of p t for the 
semicentral collisions Ni+Ni at 1.93 A GeV. In the upper panel of Fig. ^] we present by solid 
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lines the calculations employing the BR parametrization (open circles — full calculation, 
open pentagons — without V field) and by dotted lines — the calculations employing the 
KL parametrization (open squares — full calculation, open diamonds — without V field). 
One can see, that the V\ coefficient of K +, s is very sensitive to the choice of the kaon 
potential. As in the case of (p x ) vs (c.f. Fig. II the BR parametrization provides the 
best description of the data. The KL parametrization gives a too small value of v i indicating 
not enough repulsion. Neglecting the V field results in too large values of v 1 for the both 
parametrizations. 

In the lower panel of Fig. Elwe explore the relative importance of the kaon-nucleon (KN) 
scattering and of the kaon potential for the description of the directed flow. To this aim, 
we have performed three additional calculations: (i) keeping the BR kaon mean field, but 
without the KN scattering (solid line with open down triangles), (ii) without kaon mean 
field, but with the KN scattering (dashed line with open up triangles), and (iii) without 
kaon mean field and without the KN scattering (dash-dotted line with stars). We see that 
without the kaon mean field it is impossible to reproduce the measured positive v\ at small 
transverse momenta. The v% vs pt dependence for kaons is similar to the one for protons 
(c.f. Fig. IT2*j) in this case. The kaon mean field alone gives already the correct value of V\. 
The KN scattering reduces v\ slightly. This is expected, since the KN collisions should make 
kaons to "flow" together with nucleons. 

Figs. El IT5l and IT6l show the azimuthal distributions of K +7 s at midrapidity from semi- 
central collisions Ni+Ni at 1.93 A GeV, Au+Au at 1 A GeV and Au+Au at 1.5 A GeV, re- 
speetively, in eomparison to the data tan Ref, QQ. In order to q nant lty an anisotap. 
of the azimuthal distributions we have performed a fit of these distributions as in Refs. 



24, 



29j : 

^(0) oc 1 + 2vi cos(0) + 2v 2 cos(20) . (12) 
The elliptic flows v 2 for the data and for different calculations are collected in the Table |H] 

E21. 

The experimental data for all three systems reveal a pronounced out-of-plane (v 2 < 0) 
emission of K +, s. Overall, the calculations with the BR parametrization of the kaon mean 
field (solid lines with open circles) provide the most reasonable description of the data on 
the azimuthal distributions. The calculations with the KL parametrization (dotted lines 
with open boxes) produce not enough anisotropy. 
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In order to understand the mechanism of the out-of-plane K + enhancement better, we 
have performed additional calculations by switching-off various effects. The calculations 
without the space component of the kaon vector field V (solid lines with open pentagons for 
the BR parametrization and dotted lines with open diamonds for the KL one) differ only 

2311 on the influence 



very slightly from the full calculations, contrary to the findings of Ref. 
of the Lorentz force on the kaon squeeze-out. We recall, that for the in-plane flow (c.f. 
Figs. (fTT |) . (JT5j) ) there is a strong influence of the V field, in agreement with the results of 
Ref. 

The squeeze-out of particles can be caused either by a dynamical focusing due to a 
repulsive mean field j3] or by a shadowing of in-plane emitted particles by spectator pieces 



41 1 . Shadowing implies a dominant role of scattering and/or absorption of the particles on 
nucleons. What determines the kaon squeeze-out: the kaon mean field or the KN scattering? 
To answer this question, as in the case of the directed flow, we present in the lower panels of 
Figs. SUmini the results (i) with the BR parametrization, but without the KN scattering, 
(ii) without kaon mean field, but with the KN scattering, and (iii) without kaon mean 
field and without the KN scattering, (the meaning of the lines is the same as in the lower 
panel of Fig. discussed above). In the calculation (iii) there is no squeeze-out signal. 
The calculations (i) and (ii) both produce the squeeze-out. However, in the first case the 
signal is stronger than in the second case. Futhermore, the calculation (i) gives almost the 
same azimuthal distribution and the elliptic flow v% as the full calculation including the BR 
parametrization and the KN scattering. In the case (i) the mechanism of the out-of-plane 
enhancement can be only the dynamical focusing, while in the case (ii) only the shadowing is 
active. Thus, in our calculations the squeeze-out of K +, s is caused mainly by the mechanism 
of the dynamical focusing by the repulsive mean field. 

For comparison with kaons we, finally, present in Fig. El the azimuthal distributions of 
the 7r + 's at midrapidity from the semicentral collisions Au+Au at 1.5 A GeV and Ni+Ni 
at 1.93 A GeV. The corresponding elliptic flows v-i — obtained in a similar fit procedure as 
discussed above in the case of kaons — are given in the Table IIHI Since a pion potential 
is not included in our calculations, we can not expect a good agreement with the data j^l . 
Indeed, in the case of Au+Au at 1.5 A GeV we underestimate the measured squeeze-out 
signal. The same result was observed earlier in Ref. (2] for the case of Au+Au at 1 A GeV. 
For the lighter system and higher energy (Ni+Ni at 1.93 A GeV) the influence of the pion 
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mean field should be smaller and here our calculation is able to reproduce the data 



29]. 



IV. SUMMARY AND DISCUSSION 

We have studied the K + production and propagation in heavy-ion collisions at the SIS 
energies on the basis of the BUU model [3C|, |3l|, |32|. The model propagates explicitly a 
large set of the N* and A resonances and includes the two main kaon production channels 
which are relevant ones at the beam energies of 1-2 A GeV : BB — > BYK and ttB — > KY. 
The potentials of a nonstrange baryon B, of a hyperon Y and (optionally) of a kaon K 
were included in the cross sections via shifts of the thresholds. For the kaon potential we 
lave used the BR and the KL parametrizations, which have been already applied in Ref. 



11 1 . The KL parametrization can be derived from the lowest-order approximation to the 
chiral Lagrangian The BR parametrization Q], in distinction to the KL one, takes into 
account the in-medium pion decay constant in the kaon vector field and also has a larger 
kaon-nucleon sigma term. 

First, we have looked at the K + yields and spectra. The secondary processes RN — > BYK 
and 7rB — > KY have a big contribution to the kaon production (see the Table H] and Fig.Hk). 
This makes the kaon yield very sensitive to the in-medium reduction of the NN <-> NR cross 
sections. In our "standard" calculation, which includes the SM nucleon mean field, the in- 
medium NN NR cross sections |22j and the kaon mean field in the BR parametrization, 
we have obtained a very good description of the kaon multiplicities for both heavy and light 
colliding systems in all the SIS energy region, in agreement with the QMD calculations 
of Ref. |3[- For a light system C+C at 2 A GeV the p/ a b-spectra are reproduced quite 
well, which is also shown in Ref. However, for heavy systems Au+Au and Ni+Ni the 

calculated pi a b~, and m^-spectra are somewhat too hard. The repulsive kaon potential 
reduces the yield of soft kaons and leaves the yield of hard kaons practically unchanged. 
Thus, the kaon potential makes the kaon spectra harder. In most cases the slopes of the 
K + spectra are reasonably well described by the calculation without kaon potential. But 
this calculation drastically overestimates the kaon multiplicity. Moreover, the slopes of the 
E%fa- spectra for the central collisions of the heaviest system Au+Au at 1.5 A GeV can be 
only described by using the repulsive kaon potential (Fig. |7j). 

Second, we have studied the collective in-plane and out-of-plane flows of A" + 's. The (p x ) 
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vs and v i vs p t dependencies for kaons emitted from semicentral Ni+Ni collisions at 
1.93 A GeV are very sensitive to the choice of the kaon mean field and clearly favour the 
BR parametrization, as was also recently shown in Ref. The azimuthal distributions 

of kaons at midrapidity are also sensitive to the choice of the kaon mean field. The KN 
scattering alone gives a too weak squeeze-out signal with respect to the data. For the 
heaviest measured system Au+Au at 1 and 1.5 A GeV, the data on the elliptic flow v-i are 
best described by using the BR parametrization. For a lighter system Ni+Ni at 1.93 A GeV 
the calculation with the KL parametrization reproduces the data on v-i best, but the BR 
one is also consistent with the data within the experimental errorbars. 

The Lorentz force caused by the space components V of the kaon vector field does not 
significantly influence the azimuthal distributions of kaons at the midrapidity. However, in 
the case of the in-plane flow the Lorentz force is found to contribute very strongly [lj] (see 
also our Figs. llipf3|) . Thus, the kaon azimuthal distributions at midrapidity probe, basically, 
the static kaon potential. 

Some comments are also in order with regard to our results of a recent benchmark test 
of the transport codes |38j]. Our results in |38j were obtained with an enforced A(1232) 
lifetime of ft/ 120 MeV and with the the vacuum NN NR cross sections. In our standard 
calculations, however, we use the A(1232) lifetime of h/T &{M), where Fa(M) is the mass 
dependent width of the A(1232)-resonance j^J, and the in-medium NN NR cross sec- 

381 has led to enhanced 



tions. In particular, using the vacuum NN <-> NR cross sections in 
pion and kaon yields in our calculations with respect to most other transport calculations. 
However, our results l38|] on ti^'s and K +, s agree with those obtained with the BUU code 
of Refs. which also propagates explicitly a large set of the baryonic resonances 

with M < 2 GeV. Thus, our present study demonstrates clearly, that the effect of the higher 
baryonic resonances on the pion and kaon production is counterbalanced by the in-medium 
reduction of the NN <-» NR cross sections, which produces a satisfactory agreement with 
experimental data. 

In conclusion, the BR parametrization of the kaon mean field provides the best overall 
description of the K + observables at SIS energies. The remaining problem lies in the some- 
what too hard p£ a &-spectra for the Au+Au collisions. These spectra, however, are sensitive 
not only to the choice of the kaon mean field, but also to the kaon production cross sections 
in nuclear medium, which are still rather ambiguous. The in-medium calculation for the 



17 



channel ttB — > YK has been performed in Ref. 62| indicating a reduction of the cross 
section at a finite baryon density. Such an in-medium calculation is still needed, however, 
for another important channel BB — > BYK, basing e.g. on the model of Ref. 43]. 
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TABLE I: Contribution of the different kaon production channels to the total kaon number from 
central (6 = fm) collisions of the various systems. 



Channel Au+Au, 0.96 A GeV Au+Au, 1.48 A GeV Ni+Ni, 1.93 A GeV 

RN — > BYK 51% 47% 36% 

NN — ► BYK 12% 25% 47% 

itB^KY 25% 17% 10% 

RR — > BYK 8% 6% 3% 

ttB^KKN 3% 3% 2% 

BB — > NNKK 1% 2% 2% 
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TABLE II: The elliptic flow V2 obtained by the fit of the K + azimuthal distributions by a function 
oc 1 + 2v\ cos(4>) + 2v2Cos(2(p). The data are from Refs. 24. 12S | . 





Au+Au, 1 A GeV 


Au+Au, 1.5 A GeV 


Ni+Ni, 1.93 A GeV 


Exp. 


-0.110 ± 0.011 


-0.08 ± 0.02 


-0.04 ± 0.02 


BR 


-0.100 ± 0.007 


-0.082 ± 0.004 


-0.059 ± 0.002 


BR w/o V 


-0.092 ± 0.005 


-0.098 ± 0.006 


-0.076 ± 0.001 


BR w/o KN scatt. 


-0.088 ± 0.011 


-0.068 ± 0.006 


-0.061 ± 0.003 


KL 


-0.055 ± 0.005 


-0.062 ± 0.003 


-0.037 ± 0.004 


KL w/o V 


-0.080 ± 0.005 


-0.062 ± 0.003 


-0.041 ± 0.002 


w/o pot. 


-0.048 ± 0.005 


-0.040 ± 0.004 


-0.020 ± 0.002 


w/o pot. w/o KN scatt. 


-0.015 ± 0.005 


-0.017 ± 0.005 


-0.002 ± 0.002 


TABLE III: The elliptic 


flow V2 obtained by the fit of the ir + azimuthal distributions by a function 


oc 1 + 2vi cos(cp) + 2v2Cos{2(j)) . The data are from Ref. |29j]. 






Au+Au, 1.5 A GeV 




Ni+Ni, 1.93 A GeV 


Exp. 


-0.15 ± 0.01 




-0.04 ± 0.01 


BUU 


-0.063 ± 0.006 




-0.041 ± 0.008 
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FIG. 1: Upper panel: the kaon potential at zero momentum as a function of the nuclear matter 
density. Lower panel: the scalar density vs the nuclear matter density. 
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FIG. 2: Number of kaons (K +, s and K 0, s) as a function of time for central (b = fm) collisions 
Au+Au at 0.96 A GeV (a), Au+Au at 1.48 A GeV (b) and Ni+Ni at 1.93 A GeV (c) for the different 
production channels: RN -> BYK — dashed lines, ttB -»■ KY — dotted lines, NN -> - 
thin solid lines and RR — ► BYK — dash-dotted lines. Total kaon number is represented by thick 
solid lines in the panels a,b,c. Panel d shows the time evolution of a central baryon density for the 
central collisions Au+Au at 0.96 A GeV — solid line, Au+Au at 1.48 A GeV — dashed line and 
Ni+Ni at 1.93 A GeV — dotted line. 
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FIG. 3: Pion (upper two lines) and kaon (lower two lines) multiplicities per projectile nucleon in 
the C+C (highest and lowest lines) and in the Au+Au (two lines in the middle) systems vs the 
projectile energy. Pion multiplicity includes all pions and is scaled by a factor of 10~ 2 . Data are 
from Ref . Ul . 
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FIG. 4: Inclusive K + production cross section as a function of the laboratory momentum for 
Au+Au collisions at 0.96 A GeV at various laboratory angles. Errorbars on the calculated results 
are statistical. See text for further details of calculations. Data are from Ref. 41. 
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FIG. 5: Inclusive K + production cross section as a function of the laboratory momentum for C+C 
at 2 A GeV at 6i ab = 32° and 40°. Data are taken from Ref. U. 
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FIG. 6: Differential K + production cross section at @i a b = 40° from Au+Au collisions at 1.5 A 
GeV for different centrality bins as a function of the kinetic energy in the cm. system. Data points 
are taken from Ref. 28] and correspond to the following centrality bins (from top to bottom with 
decreasing centrality): 5%, 15%, 15%, 25% and 40% of the reaction cross section. Corresponding 
calculations are performed within the following impact parameter regions (from the upper to lower 
lines): 6 < 3 fm, 6 = 4 — 6 fm, 6 = 7 — 8 fm, 6 = 9 — 10 fm and 6 = 11 — 14 fm. Spectra are scaled 
by the factors 10 5 , 10 3 , 10 2 , 10 1 , 10° from top to bottom. 
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FIG. 7: Upper panel: inverse slope parameter of the K + c.m. kinetic energy spectrum as a function 
of the number of participating nucleons A par t for the system Au+Au at 1.5 A GeV. Lower panel: 
K + multiplicity per participating nucleon vs A par t for Au+Au at 1.5 A GeV. The data are from 



Ref. 



28 
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FIG. 8: Transverse mass spectra of K +, s from the Ni+Ni collisions at 1.93 A GeV for central 
events (left panel) and peripheral events (right panel) in comparison to the data from Ref. 



29| . The 

spectra are extracted in the following cm. rapidity intervals (from top to bottom): —0.69 < y cm < 
-0.54, -0.54 < y cm < -0.39, -0.39 < y cm < -0.24, -0.24 < y cm < -0.09, -0.09 < y cm < 0.06. 
The scaling factors 10 4 , 10 3 , 10 2 , 10 1 and 10° are applied to the spectra from top to bottom. 
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FIG. 9: K + cm. rapidity distributions from the central (left panel) and p erip heral (right panel) 
collisions of Ni+Ni at 1.93 A GeV in comparison to the the data from Ref. 26j]. The distributions 
for the peripheral collisions are multiplied by a factor of 4. 
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FIG. 10: Mean transverse momentum projected on the reaction plane vs the normalized rapidity 
for A hyperons (upper panel) and for protons (lower panel) for Ni+Ni collisions at 1.93 A GeV 
with b < 4 fm. The particles are selected in the transverse momentum ranges pt/m\ > 0.5 for A's 
and pt/mp > 0.5 for protons. The Central Drift Chamber (CDC) angular cuts (30° < 6j a & < 150°) 
are taken into account in calculations. The data are from Ref. 12511. 
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FIG. 11: K + mean transverse momentum projected on the reaction plane vs the normalized 
rapidity for Ni+Ni at 1.93 A GeV, b < 4 fm compared with the data from Ref. |l8]. Only high 
transverse momentum kaons (pt/rriK > 0.5) were selected according to |la |. The CDC angular cuts 
(39° < Qiab < 150°) and the upper limit of the laboratory momentum (0.5 GeV/c) up to which a 
K + can be identified |l8| were also taken into account in calculations. The data points and the 
calculated curves at 

y(o) > o are obtained by reflection from those at < 0. 
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FIG. 12: Proton directed flow vs transverse momentum for collisions Ni+Ni at 1.93 A GeV with 
b < 2.5 fm. Protons are selected in the rapidity range -1.2 < y (0) < -0.65. Data are from Ref. 
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FIG. 13: K + directed flow vs transverse momentum for collisions Ni+Ni at 1.93 A GeV with 
b < 2.5 fm. Kaons are selected in the rapidity range -1.2 < < -0.65. Data are from Ref. 
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FIG. 14: K + azimuthal distributions for semicentral (6 = 4 — 6.5 fm) Ni+Ni collisions at 1.93 A 
GeV. Kaons are selected in the rapidity range \Y^ \ < 0.4 and in the transverse momentum range 
Pt = 0.2 — 0.8 GeV/c. Different calculations are explained in the text. Data are from Ref. j^ . 
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FIG. 15: K + azimuthal distributions for semicentral (6 = 5 — 10 fm) Au+Au collisions at 1 A 
GeV. Kaons are selected in the rapidity range \ Y^\ < 0.6 and in the transverse momentum range 
Pt = 0.2 — 0.8 GeV/c. Different calculations are explained in the text. Data are from Ref. j^ . 
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FIG. 16: K + azimuthal distributions for semicentral (b = 6 — 10 fm) Au+Au collisions at 1.5 A 
GeV. Kaons are selected in the rapidity range \Y^ \ < 0.4 and in the transverse momentum range 
Pt = 0.2 — 0.8 GeV/c. Different calculations are explained in the text. Data are from Ref. 
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FIG. 17: Azimuthal distributions of 7r + -mesons produced in semicentral (b = 6 — 10 fm) Au+Au 
collisions at 1.5 A GeV (upper panel) and semicentral (6 = 4 — 6.5 fm) Ni+Ni collisions at 1.93 
A GeV (lower panel). Pions are selected in the rapidity ran ge jyCQ) | < 0.4 and in the tr ansverse 
momentum range pt = 0.2 — 0.8 GeV/c. Data are from Ref. 29]. 



